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profile  with  a  Laser  Doppler  Velocimeter  in  a  fully  developed 
turbulent  pipe  flow  of  a  polymer  drag  reducing  solution.  A 
one-equation  turbulence  model  has  been  developed  by  using  a 
damping  factor  similar  to  that  of  Mizushima  and  Usui. 
Predictions  by  the  model  agree  well  with  experimental  data 
for  the  mean  velocity  profile  and  for  the  increase  of  the 
length  scale  in  polymer  solutions.  Predictions  are  also 
compatible  with  the  streamwise  turbulent  intensity  measure¬ 
ments  and  show  similar  trends  to  the  measured  intensity 
distributions  in  Newtonian  fluids  and  polymer  solutions. 

In  highly  reduced  drag  region,  predictions,  together  with 
the  measurements,  indicate  considerable  effect  of  Reynolds 
number  below  Re  =  4xl04  on  the  turbulent  intensity  profile. 
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I .  INTRODUCTION 

Considerable  progress  has  been  made  in  recent  years  toward  the 
understanding  of  wall-bounded  turbulent  shear  flows  in  Newtonian 
fluids,1  whereas  the  same  cannot  be  said  for  the  same  type  of  polymer 
flows  which  are  of  considerable  interest  because  of  the  drag  reduction 
phenomenon.  Some  measurements  on  the  mean  and  fluctuating  velocity 
profiles  and  of  the  large  scale  structure  in  turbulent  polymer  flows 
are  available  and  are  discussed  in  a  recent  review  by  Berman.2 
Measured  mean  velocity  profiles  among  investigators3 ’ 4  * 5  agree  in 
several  aspects  and  general  features  of  the  profile  have  begun  to 
emerge.  In  the  laminar  sublayer  region,  the  velocity  profile  obeys 
u+  =  y+,  as  a  Newtonian  fluid.  A  change  in  the  velocity  profile  is 
observed  in  the  buffer  layer  which  extends  into  the  logarithmic 
region.  At  the  maximum  drag  reduction  point,  the  logarithmic  region 
disappears.  Concerning  the  turbulent  intensity  profiles,  the  few 
available  experimental  data4’ 5,6  all  conflict  on  the  magnitude  of 
the  peak  intensity.  Rudd,6  using  a  pipe  with  a  square  cross  section, 
observed  a  large  increase  in  the  axial  intensity  distribution  for  a 
polymer  solution  whereas  no  other  experimental  measurements  of  polymer 
flows  in  channels  and  circular  pipes  show  such  a  large  increase. 
Measurements  by  Mizushima  and  Usui5  and  Reischman  and  Tiederman4 
are  closer  to  each  other  than  to  the  data  of  Rudd6 .  Yet  they 
disagree  on  whether  the  peak  intensity  is  higher  in  a  polymer 
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solution  than  in  a  Newtonian  fluid.  In  this  paper,  a  one-equation 
turbulence  model  is  developed  and,  together  with  some  experimental 
measurements,  an  attempt  is  made  to  study  the  general  features  of 
turbulence  in  polymer  solutions.  An  explanation  for  the  conflict 
between  the  two  turbulence  measurements4’5  in  polymer  solutions  is 
offered. 

II.  EXPERIMENTS 

The  pipe  flow  facility  used  in  the  present  experiment  and  a 
detailed  view  of  the  test  section  are  shown  in  Figure  1.  The  straight 
inlet  pipe  length  is  approximately  160-pipe  diameters.  The  pressure 
drop  was  measured  by  a  ultra  low  range  differential  pressure  transducer 
over  a  75-inch  length  across  the  test  section.  Linearity  of  the 
pressure  drop  in  the  inlet  section  of  the  pipe  was  checked  for  both 
solvent  and  polymer  solutions.  For  velocity  measurements,  a  Laser 
Doppler  Velocimeter  was  employed.  The  measuring  section  of  the  pipe 
is  constructed  with  a  thin  polyester  film  (51  ym  thick)  similar 
to  that  of  Mizushima  and  Usui.5  This  construction  enables  velocity 
to  be  measured  close  to  the  pipe  wall  and  the  smallest  dimension  of 
the  measuring  volume  to  be  aligned  in  the  direction  of  the  steepest 
velocity  gradient,  minimizing  the  velocity  biasing7  and  the  broadening 
of  Doppler  signals  due  to  the  velocity  gradient.8  With  the  analysis 
of  laminar  flow  measurements,  the  measuring  volume  was  found  to  be 
about  0.05  x  0.05  x  0.250  (mm3),  and  the  error  due  to  the  finite 
transit  time  broadening  was  less  than  2.2  percent. 
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The  polymer  used  was  polyethylenoxide  (Union  Carbide  WSR-301)  with 
an  intrinsic  viscosity  measured  to  be  12.5  dl/g.  The  flow  rate  was 
measured  with  a  precision  flow  meter  which  was  calibrated  for  both  water 
and  polymer  solutions  prior  to  the  experiment  by  using  a  weighing  scale 
and  a  stop  watch.  At  present,  due  to  the  capacity  of  this  meter,  the 
highest  Reynolds  number,  based  on  the  pipe  diameter,  achievable  for  the 
experiments  was  limited  to  around  12,000. 


III.  TURBULENCE  MODEL  AND  THE  GOVERNING  EQUATIONS 

Because  of  the  success  of  the  algebraic  turbulence  model  of 
Mizushima  and  Usui5  applied  to  flows  with  drag  reduction,  a  simple 
one-equation  turbulence  model  which  prescribes  the  length  scale  and 
considers  the  transport  of  turbulent  kinetic  energy  was  developed. 

The  model  is  a  modification  of  the  Hassid  and  Poreh9  scheme  by  adopting 
the  Mizushima  and  Usui5  damping  factor  and  by  modifying  the  dissipation 
term  to  account  the  polymer  effects  in  the  buffer  zone.  The  eddy 
diffusivity  is  expressed  by 


' 

' 

1  -  exp 

-  Rt  |-a  +  (o^+l^ 

% 

, 

- 1 

J 

(1) 


The  damping  factor  is  similar  to  the  form  used  by  Mizushima  and  Usui , 5 
except  that  the  Van  Driest  constant,  A+,  and  y+  are  replaced  by  the 
constant,  A^,  and  the  Reynolds  number  of  turbulence,  Rt: 


k%  L 
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since  the  Van  Driest  mixing  length  model  is  found  to  give  unsatisfactory 
results  for  low  Reynolds  number  flows.  The  value  of  a  is  determined  in 
exactly  the  same  way  as  Mizushima  and  Usui,5  namely, 

2 


•  ' 

2A 

'U* 

V 

26 

4  - 

(3) 


where  A,  V,  U*  are  the  relaxation  time  for  the  Maxwell  model,  the 
kinematic  viscosity  and  the  friction  velocity,  respectively.  The 
constant  A  is  determined  from  the  relaxation  time  for  simple  laminar 
shear  flow,  A^,  as 


\  1  •  34 

=  3.76  x  108  W  ,  (4) 

Kl 


where  W  is  the  Weisenberg  number  *  A^  U/2R,  U  is  the  mean  velocity  and 
R  is  radius  of  the  tube.  The  calculation  of  A^  was  performed  using 
the  equation  derived  by  James  and  Acosta.10  All  coefficients  used  to 
determine  a  had  the  same  values  as  those  used  by  Mizushima  and  Usui. 5 

The  length  scale,  L,  is  prescribed  by  using  Nikuradse's  relation: 


Then  the  turbulent  energy  is  obtained  by  solving  the  transport  equation: 


Dk 

Dt 
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together  with  the  following  assumptions: 


-  u, 
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(8) 


In  the  region  close  to  the  wall,  viscous  diffusion  is  balanced  by 
viscous  dissipation  and  the  turbulent  kinetic  energy  varies  as  y2 .  The 
first  term  in  Equation  (8)  is  introduced  to  satisfy  these  two  conditions 
simultaneously,  since  the  second  term  alone  cannot  satisfy  them.  Jones 
and  Launder11  used  2v(3k2/3y)2  instead  of  2vk/L2 .  The  results  should 
not  differ  much  because  this  term  is  only  important  near  the  wall  and 
the  distribution  of  k2  is  nearly  linear  in  that  region.  Physically, 
this  term  is  interpreted  as  the  dissipation  due  to  the  large  eddies 
which  partially  submerge  in  the  viscous  sublayer  (Townsend).12  Lumley's13 
explanation  is  that  the  eddy  being  parasitic  does  not  contribute  to  the 
Reynolds  stress  and  that  its  size  and  intensity  are  determined  by  the 
mean  velocity  profile.  Since  its  scale  is  much  larger  than  the  polymer 
molecules  and  it  loses  energy  mainly  in  the  viscous  sublayer  region 
where  polymer  molecules  are  not  stretched,  the  dissipation  due  to  the 
large  eddy  will  not  be  affected  by  polymer  concentration. 

The  second  term  of  Equation  (8)  is  considered  to  be  the  dissipation 
due  to  small  eddies  through  the  energy  "cascade"  process.  At  high 
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turbulent  Reynolds  numbers,  dissipation  is  mainly  due  to  the  small 
eddies  and  large  eddy  dissipation  is  almost  negligible.  At  the 
buffer  zone,  however,  the  turbulent  Reynolds  number  is  not  large  so 
that  turbulent  energy  is  dissipated  even  in  the  smaller  wave  number 
region.  Therefore,  for  Newtonian  fluids  (a  =  0) ,  a  damping  factor 
{1  -  exp  (-A^Rt)},  is  multiplied  to  the  usual  CDk3^2/L  term.  In 
polymer  solutions,  polymer  molecules  are  stretched  in  the  buffer 
zone  and  the  stretched  polymer  molecule  size  is  comparable  to  the 
Kolmogorov  microscale.  Justification  for  the  term  a  +  (a2  +  l)^j 
in  the  damping  factor  of  the  dissipation  term  is  that,  due  to  the 
polymer,  the  proportion  of  the  two  types  of  dissipation  will  change. 

The  model  constants  C  and  CL  are  the  same  as  those  of  Wolf stein1 4 

\1  D 

who  determined  them  from  the  conditions  in  the  vicinity  of  the  wall 
where  convection  and  diffusion  of  energy  are  negligible.  The  value 
of  a  is  obtained  from  Equations  (3)  and  (4) ,  exactly  in  the  same  way 
as  Mizushima  and  Usui . 5  The  sets  of  constants  used  here  are  the  same 
as  those  of  Hassid  and  Poreh,9  summarized  as  follows: 

C  =  0.416,  CL  =  0.22, 

U  D 

A  =  0.012,  O  =  1.0  .  (9) 

M  K- 

The  turbulent-energy  equation  model  for  polymer  solutions  developed 
here  was  applied  to  fully  developed  turbulent  pipe  flow.  In  dimension¬ 
less,  universal  variables,  the  governing  equations  are: 


dU+  y+ 

a  +  e+>  ^  =  a  -  $  . 


(10) 
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(13) 

Numerical  solutions  of  the  resulting  equations  were  obtained  by  using  a 
second-order  accurate  finite  difference  scheme,  solved  with  Gauss-Seidel 
relaxation.  Boundary  conditions  for  the  turbulent  kinetic  energy  are: 

k  =  0  at  y  =  0  (14a) 

k  =  kg  +  k^r2  near  centerline  of  pipe  (14b) 

where  k^  is  a  constant  obtained  from  Equations  (13)  and  (14b) . 

In  order  to  place  many  grid  points  near  the  wall,  where  the  gradients 
are  large,  the  following  transformation  is  used: 

H  =  In  (1  +  y+) 


(15) 
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IV.  COMPARISON  OF  PREDICTION  WITH  EXPERIMENTAL  RESULTS 

Mean  velocity  profiles  are  shown  in  Figure  2.  For  the  maximum  drag 
reduction  asymptote  calculation,  the  value  of  a  was  set  equal  to  40 
following  Reference  5.  Velocity  profiles  calculated  with  the  same  a 
change  for  different  Reynolds  numbers.  The  prediction  for  the  maximum 
drag  reduction  at  a  Reynolds  number  of  14,400  agrees  well  with  the 
experimental  data  of  a  300  ppm  polymer  solution  by  Mizushima  and  Usuis 
at  the  same  Reynolds  number.  For  the  three  other  comparisons  made  for 
polymer  solutions  the  agreement  is  also  good.  For  each  of  these  cases, 
the  value  of  a  was  determined  by  using  Equations  (3)  and  (4) ,  together 
with  U.^,  U,  and  T  obtained  from  the  experiments.  All  three  predicted 
velocity  profiles  deviate  from  U+  =  y+  line  at  around  y+  of  10.  Measured 
mean  velocity  profiles  in  the  near -wall  region  (y+  <  25)  obtained  by 
Reischman  and  Tiederman1*  show  the  same  kind  of  deviation. 

The  coefficient  of  resistance  versus  Reynolds  number  is  shown  in 
Figure  3.  The  a  =  0  predicted  curve  coincides  with  Prandtl's  curve 
except  in  the  very  low  Reynolds  number  range  and,  as  can  be  seen,  the 
present  experimental  data  are  in  excellent  agreement  with  the  calculations . 
Various  calculations  made  for  different  values  of  a  are  also  shown.  The 
maximum  asymptote  predicted  with  a  =  40  differs  slightly  and  has  a 
slightly  flatter  profile  compared  with  that  of  Mizushima  and  Usui.5 
Points  A,  B,  and  C  in  Figure  3  were  obtained  from  the  experiments 
corresponding  to  the  velocity  profiles  shown  in  Figure  2.  Values  of  a 
shown  in  Figure  2  for  each  predicted  velocity  profile  agree  well  with 
the  values  which  can  be  deduced  from  Figure  3. 
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Turbulent  energy  distributions  are  shown  in  Figure  4.  The 
predictions  by  Hassid  and  Poreh9  show  a  large  increase  in  the  energy 
for  flows  with  drag  reduction  with  a  maximum  value  close  to  that 
measured  by  Rudd.6  Rudd's6  data  were  obtained  in  a  one-half  inch 
square  duct  so  that  it  is  highly  likely  that  the  data  were  contaminated 
with  the  secondary  flow  effects  of  the  corners.  This  is  seen  from 
Rudd's6  measurements  for  water  which  indicate  significantly  higher 
turbulent  intensity  than  other  experimental  results  (Reischman  and 
Tiederman,4  Eckelmann  and  Reichardt,15  Bakewell,16  and  Pennel,  et  al.).17 
The  present  measurements  for  water  are  in  closer  agreement  to  the  latter 
investigators'  results  as  shown  in  Figure  5.  Present  predictions 
simulate  the  tendency  of  the  data  obtained  by  Reischman  and  Tiederman4 
which  show  a  slight  increase  in  the  (u+)  peak  for  a  polymer  solution 
with  the  peak  occurring  at  about  y+  =  40.  The  data  of  Mizushima  and 
Usui,5  however,  show  a  decrease  in  the  (u+)  peak.  As  will  be  explained 
later,  this  may  be  attributed  to  the  lower  Reynolds  numbers  at  which 
the  experiments  by  Mizushima  and  Usui5  were  performed. 

Both  experimental  data  (References  4  and  5)  and  present  predictions 
indicate  a  decline  in  the  turbulent  intensity  in  the  region  y+  <  20. 
Figure  5  shows  the  turbulent  intensity  distribution  for  different 
Reynolds  numbers.  Calculated  results  indicate  the  Reynolds  number 
effect  is  significant  for  a  polymer  solution,  whereas  the  effect  for 
the  solvent  is  minimal.  The  experimental  data  support  this  tendency. 

The  turbulent  intensity  with  a  Reynolds  number  of  11,000  is  significantly 
higher  than  at  a  Reynolds  number  of  7,200.  The  data  points  (D  and  E) 
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correspond  to  points  D  and  E  in  Figure  3  and  are  for  high  drag  reduction 
conditions.  Predicted  Reynolds  number  effect  on  the  peak  turbulent 
intensity  and  its  location  are  shown  in  Figure  6,  together  with  the 
experimental  data  of  Figure  5.  For  Newtonian  fluids,  the  effect  is 
negligible  for  a  Reynolds  number  higher  than  1  x  104 .  For  flows  with 
large  drag  reduction,  the  effect  on  the  peak  turbulent  intensity  is 
significant  for  a  Reynolds  number  up  to  around  5  x  104  and  there  is  a 
Reynolds  number  above  which  the  turbulent  intensity  is  higher  than  that 
of  the  solvent  and  below  which  the  intensity  is  lower.  Experimental 
data  taken  at  Reynolds  numbers  of  7,200  and  11,000  follow  the  same 
trend  as  the  predictions  indicate.  On  the  location  of  the  peak 
turbulent  intensity,  on  the  other  hand,  the  Reynolds  number  effect  is 
small.  The  difference  of  the  (u+)  peak  for  polymer  solutions  between 
the  data  by  Mizushima  and  Usui5  and  that  by  Reischman  and  Tiederman4 
may  be  attributed  to  the  Reynolds  number  effect,  since  the  former  is 
taken  at  a  Reynolds  number  range  between  9,500  and  35,000  and  the 
latter  between  20,000  and  52,000. 

The  predicted  increase  in  the  length  scale  for  a  polymer  solution 
is  compared  with  experimental  data  in  Figure  7.  Length  scale  measured 
in  the  experiments  was  the  mean  spacing  between  the  low  speed  streaks. 

A  good  measure  of  this  spacing  should  be  the  thickness  of  the  laminar 
sublayer.  Therefore,  a  laminar  sublayer  thickness  based  on  the  maximum 
of  the  turbulent  energy  production  is  used  for  the  prediction.  The 
predicted  length  scale  increase  agrees  well  with  the  data  by  Donahue, 
Tiederman,  and  Reischman,18  and  by  Oldaker  and  Tiederman.19 
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Figures  8  and  9  show  the  turbulent  energy  balance  predicted  by  the 
model  for  a  Newtonian  fluid  and  for  a  polymer  solution  in  the  maximum 
drag  reduction  region  (a  «  40).  A  significant  contrast  between  the  two 
is  that,  for  a  flow  with  maximum  drag  reduction,  the  dissipation  due  to 
the  second  term  in  Equation  (8)  is  dominant  except  in  the  region  close 
to  the  wall,  i.e.,  y+  <  10,  and  total  dissipation  is  mainly  due  to  the 
second  term  which  is  defined  as  an  energy  cascade  type  of  dissipation. 
For  a  Newtonian  fluid,  the  first  term  in  Equation  (8),  i.e.,  the 
dissipation  due  to  the  large  eddies,  plays  a  significant  role  in  the 
total  dissipation,  whereas  the  second  term  does  not  dominate  the 
dissipation  until  y+  becomes  greater  than  around  30. 

Although  the  model  adopted  is  extremely  simple,  a  rough  qualitative 
argument  may  be  made  on  the  spectrum  of  turbulence  from  Figures  8  and  9. 
The  dissipation  associated  with  the  mean  velocity  gradient  is  expected 
to  be  at  relatively  low  frequency,  while  the  cascade  type  of  dissipation 
occurs  at  a  relatively  high  frequency.  For  a  Newtonian  fluid,  Figure  8 
shows  that  the  higher  frequency  content  will  increase  with  the  distance 
from  the  wall  up  to  y+  ~  30.  The  experimental  data  by  Bakewell,15  which 
shows  the  spectra  of  the  streamwise  velocity  fluctuations  at  y+  =  10, 

20,  and  40,  confirm  this  trend  (see  Figure  10).  For  polymer  solutions, 
there  is  no  spectral  data  available.  The  author  has  attempted  to 
obtain  these  data  but  the  small  difference  in  the  change  of  frequency 
content,  as  seen  in  Figure  10,  and  the  spectrum  error  inherent  to  the 
Laser  Doppler  Velocimeter  prevented  any  definite  conclusion  from  being 
formed . 
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Calculation  results  indicate  the  cascade  type  of  dissipation 
dominates  polymer  flows  in  the  maximum  drag  reduction  region.  It 
should  be  noted,  however,  this  does  not  imply  that  fluctuating 
velocities  in  polymer  flows  have  a  higher-frequency  content  than 
those  for  a  Newtonian  fluid.  It  simply  means  that  in  polymer  flows 
the  energy  is  cascaded  to  small  eddies  before  it  dissipates.  In  fact, 
the  spectral  analysis  of  the  fluctuating  velocity  shows  the  Newtonian 
flow  contains  higher-frequency  fluctuation  than  polymer  solutions. 

The  prediction  that  the  dissipation  associated  with  the  cascade 
dominates  in  the  maximum  drag  reduction  region  contradicts  the  argument 
by  Lumley13  who  predicted  the  turbulence  to  consist  primarily  of  large 
eddies  in  that  region.  More  experimental  investigations  are  needed  to 
shed  light  on  this  point. 


V.  CONCLUSION 

1.  Laser  Doppler  Velocimeter  measurements  were  made  of  the  mean 
and  fluctuating  velocity  profile  in  a  fully  developed  turbulent  pipe 
flow  of  polymer  drag  reducing  solutions. 

2.  A  one-equation  turbulence  model  has  been  developed,  as  an 
extenstion  of  the  Mizushima  and  Usui5  algebraic  model.  The  dissipation 
term  in  the  kinetic  energy  equation  is  divided  into  two  terms,  i.e., 
one  associated  with  the  changes  in  wall  shear  and  the  other  with  the 
energy  cascade.  The  polymer  effects  are  accounted  for  by  changing  the 
proportion  of  the  two  types  of  dissipation. 
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3.  Predictions  by  the  model  agree  well  with  experimental  data 
for  both  mean  velocity  and  turbulence  intensity  profiles,  as  well  as 
for  the  increase  of  the  length  scale  in  the  polymer  solutions. 

4.  In  the  highly  reduced  drag  region,  the  effect  of  Reynolds 
number  below  Re  =  4  x  104  on  the  turbulence  structure  seems  important. 
Conflicting  results  between  Reischman  and  Tiederman4  and  Mizushima  and 
Usui5  over  the  axial  turbulent  intensity  profile  may  be  attributed  to 
this  Reynolds  number  effect. 

5.  More  experimental  investigations,  including  the  study  of 
Reynolds  number  effects,  are  needed  to  clarify  the  turbulence  structure 
for  flows  with  drag  reduction. 
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Comparison  of  Calculated  and  Measured  Mean  Velocity 


Turbulent  Energy  Distribution 


Figure  5  Turbulent  Intensity  Distribution 
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PERCENTAGE  DRAG  REDUCTION 

Figure  7.  Calculated  and  Measured  Changes  in  Length  Scale 


Figure  8.  Turbulent  Energy  Balance  (Newtonian  Fluid) 


POLYMER  SOLUTION 


Turbulent  Energy  Balance  (Polymer  Solution) 
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Figure  10.  Spectra  of  the  Streamwise  Fluctuating  Velocities 
Outside  the  Viscous  Sublayer 
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